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Abstract. In this paper, we study the global subsonic irrotational flows in a multi-dimensional 
(?i > 2) infinitely long nozzle with variable cross sections. The flow is described by the inviscid 
potential equation, which is a second order quasilinear elliptic equation when the flow is subsonic. 
First, we prove the existence of the global uniformly subsonic flow in a general infinitely long 
nozzle for arbitrary dimension for sufficiently small incoming mass flux and obtain the uniqueness 
of the global uniformly subsonic flow. Furthermore, we show that there exists a critical value of 
the incoming mass flux such that a global uniformly subsonic flow exists uniquely, provided that 
the incoming mass flux is less than the critical value. This gives a positive answer to the problem 
of Bers on global subsonic irrotational flows in infinitely long nozzles for arbitrary dimension 
[3. Finally, under suitable asymptotic assumptions of the nozzle, we obtain the asymptotic 
behavior of the subsonic fiow in far fields by a blow-up argument. The main ingredients of our 
analysis are methods of calculus of variations, the Moser iteration techniques for the potential 
equation and a blow-up argument for infinitely long nozzles. 



1. Introduction 

This paper is devoted to the existence and the uniqueness of global subsonic flows for the Eu- 
ler equations for steady irrotational compressible fluids. Our focus is on the global subsonic flows 
in a general multi-dimensional infinite nozzle, which is an important subject in gas dynamics 
{see mm l8j[12Ji22j). 

Consider the steady isentropic compressible Euler equations 

div( pu) = 0, in J7, 

(1.1) 

div{pu ® u) + Vp = 0, in Q,, 

where p, u = (ni,...,u„), p represent the density, velocity, and the pressure of the fluid, 
respectively. Moreover, the pressure p = p{p) is a smooth function of p and p'{p) > 0, p"{p) > 
for p > 0. 

It is easy to derive the following so-call Bernoulli's law [5] 

u-v(^\u\^ + h{p)] = 0, (1.2) 
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fP p'{s) 

where h(p) is the enthalpy, defined by hip) = / ds. The relation (II. 2p implies that 

1 ^ 
the quantity -B(/0, = 2 1^1^ ~^ ^{p)i named Bernoulli's function, remains constant along the 

stream line in a steady isentropic flow. 

If, in addition, the flow is assumed to be irrotational, ie. the vorticity of the flow velocity 

V X ti = 0, in $7, 

then there exists a velocity potential function 93, at least locally, such that 

u{x) = Vip{x). 

In this case, the relation ()1.2p simplifies to the following strong version of the Bernoulli's law 

VB(p, |<^|2) = V Q|V(/.p + h{p)^ = 0. (1.3) 

This yields a density-speed relation for steady irrotational flows. Therefore, the density p can 
be determined by the speed |Vv3|, denoted by p Then the steady Euler equations (jl.ip 

are reduced to the following well-known scaler potential equation 

div(p(|V(^P)V(/j) = 0, in f]. (1.4) 

One of the most important parameters to the fluid dynamics is the Mack number, which is 
defined as a non-dimensional ratio of the fluid velocity to local sound speed, 

M 



where c{p) = \/p'{p) is the local sound speed. Mathematically, the second-order nonlinear 
equation ()1.4p is elliptic in the subsonic region, ie. M < 1 and hyperbolic in the supersonic 
region where M > 1. 

Subsonic flows are those in which the local velocity speed is smaller than sonic speed every- 
where, i.e. the Mach number of the flow is less than 1. Since the corresponding equations of 
subsonic flows possess some elliptic properties, problems related to subsonic flows are, in general, 
have extra-smoothness to those related to transonic flows or supersonic flows. There are many 
literatures in this field in the past decades. The first result is due to Frankl and Keldysh [15) . 
They studied the subsonic flows around a 2D finite body (or airfoil) and proved the existence 
and the uniqueness for small data by the method of successive approximations. Later on, Bers 
[T][2] proved the existence of subsonic flows with arbitrarily high local subsonic speed for the 
Chaplygin gas (minimal surface). By a variational method, Shiffman |25| j26| proved that, if 
the inflnite free stream flow speed Uoo is less than some critical speed, there exists a unique 
subsonic potential flow around a given proflle with finite energy. Shortly afterwards, Bers [5] 
improved the uniqueness results of Shiffman. Finn and Gilbarg [13] proved the uniqueness of 
the 2D potential subsonic flow about a bounded obstacle with given circulation and velocity at 
inflnity. All above the results are related to two dimensional problems. For three (or higher) 
dimensional case, Finn and Gilbarg [14) proved the existence, uniqueness and the asymptotic 
behavior with implicit restriction on Mach number M. Payne and Weinberger [23] improved 
their results soon after. Later, Dong [9] extended the results of Finn and Gilbarg [14) to any 
Mach number M < 1 and to arbitrary dimensions. Furthermore, in [10], Dong and Ou extended 
the results of Shiffman to higher dimensions by the direct method of calculus of variations and 
the standard Hilbert space method. 

All results as above (including [16)-[20)) are related to the subsonic flows past a proflle. 
Another important problem is the study of subsonic flows is the theory of global subsonic flow 
in a variable nozzles as formulated by Bers in [5]: 
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Problem 1. Find ip such that, 

dw(^p{\Vip\'^)\/Lp 

1^ = 0, 

on 

J So ol 
|V(^| < c(/)), 



0, 



mo > 0, 



in 0, 
on dQ, 

in Q, 



(1.5) 



where il. C M" is an infinitely long nozzle, niQ > is the mass flux passing through the nozzle, 
Sq is an arbitrary cross section of the nozzle, n and I are the unit outer normal of the domain 
il. and Sq, respectively (Please see Fig. [ip. 




Figure 1. Subsonic flow in a nozzle 



In the famous survey [5], Bers claimed without proof the unique solvabihty of sufficiently 
slow subsonic irrotational flows in two dimensional channel. The rigorous proof of this fact 
was achieved mathematically recently by Xie and Xin [27]. They established a very complete, 
satisfactory and systematic theory for the two dimensional subsonic flows in an infinitely long 
nozzle for potential flows, which not only solves the Problem 1 in this case, but also yields 
the existence of subsonic-sonic flows in the nozzle as limits of subsonic flows. One of the key 
ideas in [57], is to use the stream function to formulate the problem to a quasilinear elliptic 
problem with Dirichlet boundary conditions. The benefit of the stream function formation of 
the problem is that, the stream function has a priori L°° bound, and the flow region of two 
dimensional nozzle, though infinitely long, has finite "width". So one can obtain the boundary 
L°° estimate of the gradient of the stream function, Vip, by constructing proper barrier functions 
and the standard comparison principle for subsolution to second order elliptic equation. Similar 
approach has been applied in 3D axis-asymmetric nozzles by Xie and Xin in [28]. Furthermore, 
these ideas are also useful to study the physically more important case, subsonic Euler flows, by 
Xie and Xin in [29] (see also the generalization in [H]). However, it seems difficult to apply the 
method in [271 and j28| in general multi-dimensional (n > 3) nozzles, since the stream function 
formulation can not work in this case. Thus, we have to consider a different approach from that 
in [27] to treat the subsonic problem in multi-dimension case. 

On the other hand, since the domain of an infinitely long nozzle is differentiable homeomor- 
phism to an infinitely long cylinder which is unbounded, the nozzle fiow problems are different 
to the airfoil problems in which the domains are exterior domains. The main advantage of the 
exterior domain is that it can be transformed to a bounded domain through a Kelvin-like trans- 
formation. Then the airfoil problem can be transformed (explicitly or implicitly) to a scalar 
quasilinear elliptic problem with a bounded domain. This feature of the exterior domain plays 
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an essential role in the previous airfoil results. For instance, in [21], a Hardy-type inequality in 
the exterior domain is essential. But there is no similar Hardy-type inequality for the domain of 
nozzle flows, which is the one of the main difficulties in our case. For more detailed discussions, 
we refer to [30] . 

The main purpose of this paper is to study subsonic flows in general multi-dimensional 
(n > 2) infinitely long nozzles. First, we formulate a subsonic truncated problem, which is 
a uniformly elliptic equation in a bounded domain. Moreover, we prove the existence of the 
weak solution to the truncated problem by a variational method, and use the approximated 
variational problems in bounded domains to approximate the original Problem [TJ To realize this 
procedure, some uniform estimates are needed to show that the approximated solutions converge 
to the ones of the original Problem [TJ However, one can not expect to get the uniform boundary 
gradient estimate of by the classical barrier function argument, since the potential function 
is essentially unbounded, which is another main difficulty in this paper. The key observation 
here is that, though the potential function ip is unbounded, the average of V(/j is uniformly 
bounded (see the estimate ()3.6p ). Using this fact and the uniform ellipticity, we prove the "local 
average estimate" which states that the average estimate implies the local average of the gradient 
is uniformly bounded (see (|3.23p for details). That is, is locally bounded. Then, it 
is easy to get the bound of V(/9 by the standard Moser iteration. With this key estimate 
of uniformly L°° bound of we establish the existence of the subsonic flows in an infinitely 
long nozzle for arbitrary dimensions for suitable small incoming mass flux, including the two 
dimensional case in [27]. Next, we show that the global uniformly subsonic flow is unique. The 
proof is based on considering the linear equation satisfied be the difference of two solutions of 
the nonlinear potential equation. Moreover, we prove the existence of the critical incoming mass 
flux for subsonic flows. Finally, with the additional asymptotic assumptions on the nozzle at the 
far field, we obtain some asymptotic behaviors of the subsonic flow at the far field by a blow-up 
argument. 

Before stating the main results in this paper, we first give the following assumptions on the 
nozzle. 

Basic assumptions on 0. There exists an invertible C^'" map T : 17 ^ C : x ^ y 

satisfying 

T{d^) = dC, 

' For any kGR, T{n n {x„ = k}) = B{0, 1) x {y„ = k}, (1-6) 

where K is a uniform constant, C = B{0, 1) x (—00, 00) is a unit cylinder in M"", -6(0, 1) is unit 
ball in M"~^ centered at the origin, x„ is the longitudinal coordinate. 




Figure 2. Basic assumptions on 
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Asymptotic assumptions on 0. Suppose that the nozzle approaches to a cyhnder in the 
far fields, ie. 

Q n {xn = k} ^ S±, as /c — >■ ±00, (1-7) 
respectively, where S± are n — 1 dimensional, simply connected, C^'" domains. 




Figure 3. Asymptotic assumptions on Q 

Nondimensionalization of the quantities. It follows from Bernoulli's Law (jl.'2|) that in 

a potential flow the density is a given function of speed. Applying the fact ([S], [5]) that there 

exists a critical speed qcr such that the flow is subsonic if the speed is less than qcr, we can 

introduce the nondimensional velocity and density as 

u ^ V ^ p 
u = — , v = — , p = 

Qcr Qcr P[qir) 

With an abuse of the notation, we still denote the nondimesional quantities by u, v, p. Then it 
is easy to check that pq < 1 for q > and that the flow is subsonic provided that q < 1 oi p > 1. 
Our main results in this paper are stated as follows. 

Theorem 1.1. Suppose that the nozzle satisfies the basic assumptions lil.6\) . Then 

(i) there exists a positive number Mq depending only on $7, such that if niQ < Mq, then there 
exists a uniformly subsonic flow through the nozzle, ie., the ProblemUl has a smooth solution 
(f e C°°{n). Moreover, 

||Vv?(x)||ci.«(Q) < Cmo, 
where C > is a uniform constant independent of MQ,mQ, and (p. 

(ii) There exists a critical mass flux M^. < 1, which depends only on Vl, such that ifO < niQ < 
Mc, then there exists a unique uniformly subsonic flow through the nozzle with the following 
properties 

Q(mo) =sup|V(^| < 1, (1.8) 

and Q{mo) ranges over [0, 1) as mo varies in [0, Mc). 

(Hi) Furthermore, assume that the nozzle satisfies the asymptotic assumption |j.7D , then the 
flow approaches the uniform flows at the far fields, ie. 

Vip = {0,- ■ ■ ,q±), as Xn-^±oo, (1.9) 

respectively, with q± being constants determined uniquely by 

here \S±\ represents the measure of the domain S±, respectively. 

Remark 1.1. In the first statement of Theorem 11.11 it follows from the proof in Section 3 that 
one can derive an explicit form of Mq, which depends only on the nozzle In particular, it 
does not depend on the equation of the states. On the other hand, in the second statement of 
the Theorem I l.H we just give the existence of the critical mass flux Mc for a given infinite long 
nozzle. Clearly, Mq is a lower bound of Mc- 
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Figure 4. Asymptotic behaviors of the subsonic flows at the far fields 

Remark 1.2. In the proof of the uniqueness of uniformly subsonic fiows, it is not necessary to 
require the asymptotic assumption ()1.7p on the nozzle. It is quite different from the strategy in 
|27j for the 2D case, in which the proof of the uniqueness depends on the asymptotic behaviors 
of the uniformly subsonic flow in the far fields. However, in this paper, the uniqueness of the 
uniformly subsonic flow is obtained in arbitrary dimensional nozzle without the asymptotic 
assumption (|1.7p . 

This paper is organized as follows. In the next section, we introduce some necessary pre- 
liminaries. In Section 3, we prove the first statement of the Theorem 1.1. Our strategy for the 
existence of subsonic flows with small incoming mass flux can be divided into six steps: Step 1, 
truncate the coefficients of the potential equation to guarantee the strong ellipticity and truncate 
the unbounded nozzle to a series of bounded domains $7^, to formulate the approximated strong 
elliptic problems in bounded domains. Step 2, solve the approximate truncated problems by 
a direct variational method. Step 3, improve the regularity of the variational solutions to give 
the regularity. Step 4, prove the local average estimates to the gradient of the solutions. 
Step 5, obtain the classical C^'" estimate of the approximate solutions. Step 6, based on these 
key estimates, the existence of the subsonic solution to the nozzle problem for suitable small 
incoming mass flux is proved. The uniqueness of the uniformly subsonic solution is given in 
Section 4, while the existence of the critical value for incoming mass fiux is obtained in Section 
5. In the last section, we prove that the subsonic nozzle fiows approach to the uniform fiows at 
the far fields when the nozzle satisfies the asymptotic assumption ()1.7|) . 

In this paper, x, y always denote the variables in J7 and C respectively, denotes the function 
defined in Q and (p = ipoT~^ denotes the corresponding function defined in C. d, and V denote 
the derivatives with respect to x in fi, while d, and V denote the derivatives with respect to y 
in C. y4 ~ i? means 

^A<B< CA, 

with C a positive constant. 

2. Preliminaries 

In this section, we give some basic notations, definitions and facts to be used in this paper. 
2.1. Morrey theorem. 

Definition 2.1. Let Q be bounded region in M". 17 is said to be of A-type, if there exists a 
positive number A such that, for any x € 17 and < r < diam Q, 

\nnB{x,r)\ > Ar"". 
Now, we state the following Morrey theorem (see, for instance, [7]): 

Theorem 2.1. Assume that Q is of A-type, u S W^'P{VI), p > 1 and there exist constants 
K > {), 0<a<l such that, for any Br, 

I \Vu{x)\Pdx < KR^'-P+^P 
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holds. Then u € C"{Tl) and 



oscB„u < CKR" 



where C depends on n,a,p and A. 



2.2. Uniform Poincare inequality. 

Here, we prove a useful lemma of Poincare type inequality. Assume that S* is a bounded 
domain in M", and there is a constant C{n,p,U) such that the following classical Poincare 
inequality holds: 



Hy) \'dyj < C{n, p, U) \Vyu{y) \Pdy 

with / u{y)dy = 0, where C{n,p, U) depending only on n,p, [/, not on u. 
Ju 

Define a class Vk by 

Uk = { ^ I 3 an invertible smooth mapping T : Q ^ U, such that 
||r, r~^||c2,a <K< oo}. 
Then, the following useful uniform Poincare type inequality holds: 

Proposition 2.2. For any 1 <p < oo, there exists a constant C{n,p,U, K) depending only on 
n,p, U,K, such that, for any Q € Vk, 



[ \u{x)\Pdx <C{n,p,U,K) [ \Vu{x)\Pdx 
Jn Jn 



or 



\\u\\LP{n) < C{n,p,U, K)\\Vu\\lp(^^) 

holds, provided that / u{x)dx = 0. 
Jo, 

Proof. Set a = —— I uo T~^(z)dz and J = -—. Since 
\U\ Ju dy 

/ uoT~^{y)Jdy= / u{x)dx = 0, 
Ju Jn 

one has 

a|0| = / [a — uoT~^{y)) Jdy. 
Ju 

Then, by the classical Poincare inequality for p = 1 on U , one gets 

\a\\n\ = 

Ju ^ 

dy 



(2.1) 



< IIJ 



(noT ^(y) - a) Jdy 
u o T~^{y) — a 

< ||J||LooC(n,l,C/) / \S/y{uoT-\y))\dy 

Ju 

< ||J||L«.C(n,l,f/)||Vr-^||_ / \\/^uoT'\y)\dy 

Ju 

<C(n,l,C/)||Vr-i||L=o||J||i 

< C(n,l,f/)||Vr-i|Uoo||J||i 



V,uoT~\y)\''dyj^\U\^^p 

V^u{x)\Pdx]''\U\'^~p. 



1 


' ( 




J 




, In 
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Hence 



\a\\U\p < C{n,l,U)\\VT-^\\L^\\J\\L^ 



<C{n,l,U)\\\/T~'\\L^\\J\\L^ 



1 




J 




1 


p 


J 





l|Vx^/(x)||iP(f,)M 



(2.2) 



\Lp(Q.)- 



On the other hand, by the classical Poincare inequality for p on U , one has 



\u{x)\^ —dx 



< 



u 



u 



uoT-\y)fdy^ 

uoT^\y) - a\Pdyj " + \a\\U\v 
<C{n,p,U)( [ \Vy{uoT-\y))\Pdyy +\a\\U\-p 



u 



<C(n,p,f/)||Vr-i||Loo( / \VM^)\''jdx] +\a\\U\-^ 



< C{n,p, c/)||vr-i||Loo 



l|V.n(x)||iP(f,) + |a||[/|p 



<C{n,p, U)\\VT~^\\l^\\J\\lo 
which, together with (|2.2p shows 



\u{x)\P-dx) <C{n,p,U)\\VT~'\\L^ 



|V^n(x)||LP(n) + \a\\U\p, 



l|Vxn(x)|[ip(o) + |a||C/|^ 



< C{n,p, U) + C(n, 1, U) WVT'^U^ \\J\\loo 



\\^xU{x)\\LPi^^). 



Therefore, 

\\u{x)\\LP{n) 



n 

<\\J\\t 



\u{x)\^ —Jdx 



1 

1 \p 
\u{x)Y'—dx 



< (c{n,p,U) + C{n,l,U)y\VT- 



1+i 



||Vx'u(x)|IiP(!^), 



which implies the inequality (|2.ip . 

Theorem 2.3. (Uniform Poincare Inequality) For any o € M, I < p < oo, one has 



fix) 



f{x)dx 



□ 



< 



LPina,a + l) 



(2.3) 



Here 



/ f(.x)d 



{x = (xi, • • • , Xn) G ri|a < x„ < 6}, 

fix)dx, 



\^a,a-i 



1 + 1 I jCa,ti + l 

C is a positive constant depending only on n,p,Q, independent of f , a. 



Proof. This can be easily deduced from Proposition [27 



□ 
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2.3. Basic properties of 17. 

Lemma 2.4. Under the assumption il.6\) . for any < k < 2, if (p C^''^{il.), ip = ip o S 
C^'°'{C) and vise versa, and 

Similar equivalence holds for H'^ norms, s = 0, 1, 2. 

Proof. The proof follows from simple calculations, and is omitted. □ 

According to Lemma[231 we may abuse a bit of the notations by simply denoting ||'/'||c7'^.<»(r2) 
and ||(^||cfc.<i(c) by ||v||c'''<^ ll'^llc^.^i \W\\H^{n) and ||(^||iys(c) by \\lp\\s or \\^p\\s respectively. 

Lemma 2.5. Assume that 17 satisfies Then for any xq S dVt, there exists an invertible 

C^'"^ map ■ — > Bs^ : x ^ y satisfying the following properties 



T^o (U^o ^^) = K , T,o (U.o n on) = Bs, n {y„ = 0}, (2.4-1) 

CTijCTinix) = aijainiy) = 0, for X e dfl {i.e. yn = 0), 1 < j < n - 1, (2.4-2) 

11^.0,^-1^2,. <i^, (2.4-3) 

\aij{x)C,\, \a^,iy)^j\ ~ |^|, Vx G V y G B+ V ^ S M". (2.4-4) 



where Uxq is a neighbourhood of xq in M", B^^ is a ball centered at the origin with radius 5q, 

dy ' 

= Bsq n {y-ri > 0}, aij = TT^; ^0 o-iT-d C are positive numbers independent on xq € 517. 

(JXi 

Proof. By assumption (jl.6p . T(xo) = xq £ dC. Set 

Vx, =T~^By ^{xo)nc), Sxo=WondQ, s^, = By^{x^) n dc, 

then 

Suppose that x{yi, yn-i) = x{y') is the standard surface parameter of ^^^(and then, of Sxq), 
Nxoiy') is the unit inner normal vector on Sxq- Let 

h{y') = ei{y')-Nx,{y'), 1 < i < n, 
where ei, (1 < i < n) are the unit coordinate vectors. Then 

Nx,{y') = (Ai(2/'),...,An(2/')), ^iv') € C^'"- 

Define y = Tx^ [x) by 

ryi+Vn ry2+yn rVn-l+Vn 

Xi = Xi{y') + yl'"^ / •••/ Xi{si,...,Sn-i)dsids2---dsn-i, (2.5) 

■Jyi •'y2 Jyn-i 

for 1 <i <n. Since ||Aj(y') H^;!,. < C\\T, T~^'^(j2,a, there exists a (5o > independent of xq such 
that T~^{y) is well defined on and UTajg, T~^\\c^,a < K. Then, define Uxq = T~^{Bsq). 
Clearly Txo satisfies dS^l) and ([233). 

Denote the matrix ((Tjj(x)) by ^(x). For any ^ € M", 

\Aix)c\ < \A{xm\ < \A~\x)c\ < \A~\xm\ < 

Then ()2.4l -4) follows immediately. 

To prove 1)2.41 -2). we differentiate ()2.5p with respect to yj and note that y = {y , 0) for x G 50, 
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Since 



hence 



^-fo'.o) = (.,)-'fa',o) = (^fe',o))" ^di*"'") 

On the other hand, since 

<T™(y',0)--^(?y',0) = 0, l<j<n-l, 

(o"r!,i(y',0), ...,o-„„(y',0)) is the normal direction of S^o at x G S^j^, that is, (cr„i(y',0), 0)) 
is parallel to the inner normal Nxq{x) = (Ai(?/', 0), An(?/', 0)). Comparing with (j2.6p yields 

tTi*(y',0)a„i(y',0) = 0, l<j<n-l. 

□ 

Remark 2.1. Hypothesis ()1.6p is stronger than the C^'"-regularity hypothesis on Q. If one only 
assumes that G C^'"^, then C and (5o in lemma [23| in general, may depend on xq € 

Lemma 2.6. There exists a (5i > such that 
where T^^ and 6o are the same as in Lemma \2. 5\ 

Proof. Since ||r, T~^\\c2,a < C, \xi — X2I ~ \T{xi) —T{x2)\, there exists a constant 5 ~ (5o such 
that 

Taking — 2 y^^^*^^ Lemma. □ 



3. The existence of subsonic flow for small incoming mass flux 

There are two major obstacles to solve the Problem [TJ First, the ellipticity of the equation 
(|1.5p is not guaranteed beforehand, since there is no a priori bound for V(/5, the gradient of the 
solution to the Problem [TJ Second, the nozzle region is unbounded, and can not be transformed 
to a bounded domain by Kelvin-like transformations. In order to overcome these difficulties, 
we first truncate the coefficients of the equation in (jl.Sp to ensure the strong ellipticity, and 
then, truncate the domain to a series of bounded domains il^, with additional boundary 
conditions. Therefore, to solve the Problem [1] becomes to study a series of approximate strong 
elliptic problems in bounded domains and their uniform estimates, which ensure to pass the 
limit of the approximate solutions to the Problem [TJ 



3.1. A subsonic truncation and approximate solutions. 
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3.1.1. A subsonic truncation. By normalizing the equation if necessary [5], [27], one can assume 
that the critical sound speed of the flow is one. Thus, the density-speed relation ()1.3p . p = p{q^), 
is positive, sufficiently smooth and nonincreasing in q = \V{p\ G [0,1]. However, the potential 
equation is not uniformly elliptic as q approaches to 1. To guaranteed the uniformly ellipticity, 
we truncate the coefficients as follows. 

Define two functions and F{q'^) as follows 



Sis' 



Pis'), 

monotone and smooth, 
p{l - 5o), 



if < 1 - 2^0, 

if 1 - 2^0 < < 1 - 5o 
if > 1 - <5o, 



(3.1) 



and 



(3.2) 



where Jq > 0. Moreover, G(s^) is a smooth non-increasing functions and F{q'^) is a smooth 
increasing function. Set 

a^J(Vip) = Q{\S/ip\')S,j + 2Q'{\S/ip\')diipdjip. 
It is easy to check the following facts. 



F [q^) ~ q\ 



1 



<e(s^), e(s^) + 2e'(s2)s2 < C7(<5, 



C(5o) 



and there exist two positive constants A and A, such that 



(3.3) 



(3.4) 



where C{6q), X and A depend only on the subsonic truncation parameter Sq. Note that a solution 
of the potential equation derived from the new density-speed relation Q{q') is also a solution 
of the actual potential equation provided that jVc/^p < 1 — 26o- Therefore, in the end of this 
section, we will show that the solution of the truncated problem satisfies |V(/9p < 1 — 26o, as 
long as the incoming mass flux mo is suitable small. Consequently, the subsonic truncation can 
be removed. 

3.1.2. Domain truncation. Our strategy to deal with the unbounded domain here is to construct 
a series of truncated problems to approximate the Problem [1] with subsonic truncation. 

Let L > be sufficiently large. Define 

nL = {xen\\xn\<L}, = nr\{xn = ±L}, Sl = sius^. 

Consider the following truncated problem with ttiq > 0. 
Problem 2. Find a ip such that, 



div{e{\Vip\'')Vip) 

^ = 0, 
dn 

e(|v^n 

¥' = 0, 



0, 



2n dip 



mo 



dxr. 



on 



(3.5) 



on S 



L ■ 
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The additional boundary condition on S*^ implies the mass flux of the flow remains rriQ. 

Clearly, the truncated problem[2]is a strong quasilinear elliptic problem in a bounded domain. 
From now on, instead of the original Problem [H we consider a series of the truncated Problem 
[2] for any fixed sufficiently large L. With some uniform estimates of the approximate solutions, 
we can conclude that the solution of the truncated problem [2] converges to the original Problem 

m 

3.2. Truncated variational problem. 



In this subsection, we solve the truncated problem [2] by a variational method. Define 

HL = {veH\nL) : =0}. 

Then, Hl is a Hilbert space under H^-norm. The additional boundary condition on SJ^ is 
understood in the sense of traces. Define a functional J{^p) on as 

J(^)= / F(|VVp)dx-^ / i;dx', 

where F{(f') is defined by (j3.2p and x' = {xi,X2-, ■ ■ ■ , Xn-i)- The existence of solution to problem 
[2] is equivalent to the following variational problem: 

Problem 3. Find a minimizer (p G Hl such that 

J (if) = min Jiip). 



Theorem 3.1. Prohlem\^has a nonnegative minimizer ip G Hl. Moreover, 



where the constant C does not depend on L. 

Proof. Step 1. J{ip) is coercive on Hl- In fact, by Lemma[231 for any -0 G Hl, 



/ Ipdx' 


< c 


[ i'dy' 


< c 


/ / dni'dyndy' 


Jst 




Jb{o,i) 







(3.6) 



<C f \Vi)\dy <C [ \Vi;\dx ^^'^^ 

<C\nL\"^\\V^l^\\L2 
Therefore, applying (|3.7|) and Cauchy inequality yields 

JW= [ F{\Vi;\^)dx - ^ [ ^dx' 

>\[ iVV'pdx - C{mo, \S+\, lOiDllVV^Ilia 

>^\\Vi;\\l,-jC{mo,\St\,\nL\), 

which implies J('i/') is coercive. 

Step 2. The existence of the minimizer p G Hl- Since J{ip) is coercive in Hl, there is a 
minimizer sequence C Hl such that 

J{(Pn) — )• a = inf J{ip) > — oo. 

V'6-H'i 
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Then, 



<lj{o) + ^c{mo,\s+\,\nL\) 

1 



- ^^C{mo,\Sl\,\nL\). 

Therefore, there exists a subsequence, denoted by {ifn} converges weakly to some ip G and 

\\Vip\\l,<^C{mo,\St\,\nL\). 
By Fatou's Lemma, it is easy to check that 



/ F(|V^p)(ix < hminf /" F(| V(/?„p)dx. 



(3.8) 



On the other hand. 



{(pn - iffdx' <C I i'fn- <f>ydy' 



B{0,1) 



< c 



{(fin - (fi)dni<fn " 'f)dyndy' 



IB(0,1) J-L 

Jc 

J nr. 



< c 



Lfn- 'P'\ dx] I \Vtpn - V(^| dx 



as n — 7> oo. Then, 



/ \Lpn — fldx' — )• 0, as n ^ oo. 
Jst 



(3.9) 



Therefore, it foUows from ()3.8p and ()3.9p that 

>/('/5) < hminf J{ipn) = ce. 

n->oo 

i.e. 

J{ip) = min J(V') = a. 

Step 3. = max{(/7,0} is a nonnegative minimizer in H^. Indeed, since G if/,, i/?^ G if/,, 
and 

|V^+P<|V(^|2, F(|V^+p)<F(|Vv.|2), 



mo 



ip~^dx > 



\s- 



ipdx' . 



Hence, 

J(</'+) < J(¥'). 

Since (/9 is a minimizer, J((/3'^) = J{^), which imphes that ip'^ > is also a minimizer. 
Step 4. By direct computations, 

/ F{\V^\^)dx = J{ip) + / ^dx' < J(0) + ^ / V'dx' 
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It follows from ([33]) and ([321) that 

liv^lli. <\ [ F{\v^\')dx < c\^\nL\'-\\v^h2. 

That is 



I.e. 



m, 







where Smin denotes the minimal of \S'l\. □ 

Remark 3.1. The estimate ()3.6p is the key estimate for the existence of the classical solution to 
Problem [21 Indeed, the potential ip is essentially unbounded, one can not expect to get uniform 
bounds on ||V(/9||loo through as in the standard elliptic theory. 

Proposition 3.2. ip G Hi is a weak solution to the equations in i3. 5|) in the following sense: 



[ e(|V99p)V99 • Vi^dx - ^ [ ijdx' = 0, y i^eHL 

Jnr \Sr I Js+ 



(3.10) 



Proof. This is a standard variation problem. In fact, for any t € M, t > and any ^ S Hl, 
if + tip G Hl. Then, 



0<J{^ + tiP)-J{ip)= [ F{\Vip + tVi^f)-F{\V^\^)dx-^ [ ijdx'. (3.11) 

Mean value theorem yields that 

[ F(|V99 + tV^|2) -F(|V99p)(ix 

= [ [ F'{9\Vip + tVi;\'^ + \VipWl-e))de{\Vip + tViP\'^ -\Vip\'^)dx (3.12) 
Jn^ Jo 

= [ [ F'{\\/^\^ + e{t^\Vij\^ + 2tVip ■Vip))de{t^\V'ip\^ + 2tVi; ■Vip)dx. 
Jul Jo 

Since |F'(-)| < C, S/ip, Vip G L^{^l), substituting (fXT2]) into (|3TT]1 shows that 

< liminf-(J(99 + tV) - J{v)) 
t^o+ t 

= liminf / / F'{\Vip\^ + e{t^\V^\^ + 2tV^ ■Vip))de{2Vtp ■Vy^)dx 
t^o+ Jur Jo 



1 



ipdx' 



F'{\Vip\'^)de{2V^j ■Vip))dx - ^ / ipdx' 
n^Jo l^LlJst 



L 



(by Lebesgue's theorem) 

= / e(|Vv7p)V(/7- VV'dx- / i)dx'. 
JnL 

Therefore, for any -0 € Hi, 

e{\V^\'^)Vip-Vijdx- ^ [ ipdx' 

Pr I J St 



□ 
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3.3. regularity of the weak solution. We are now ready to improve the regularity of the 
minimizer (p. Indeed, one has 

Proposition 3.3. ip G Moreover, 

dp{y) 



To prove this, one needs the following estimates in Lemma [3T 
Lemma 3.4. (Interior Estimate) For any B2r{xq) C I^l, R < 6i 

Here 5i is the same number in the LemmaW. 



= 0. (3.13) 

J/n=0 



Proof. For any B2r{xq) C ^l, v S Hq (^Bs^^, h < -R, one has 

0=/ Q{\Vp\^)Vp-V{5_hv)dx = - I 5h{Q{\Vp\^)Vp)-Vvdx, (3.14) 

def 1 

where 5hv{x) = '^(^(^ + ^^^fc) — v{x)) is the k-th difference quotient, k = 1,2, ■ ■ ■ ,n. 
Set 

q = tVp^ + {l-t)Vp, ip^{x) = ^{x + hek), 

aij{q,t) = Q{f)Sij + 2Q'{q^)qiqj, a^- = aij{q) = / aij{q,t)dt. 

Jo 

Then, direct calculations give 

(5;,(e(|V(^|2)Vv?) = [ aij{q,t)dtdj{5hip) = aijdj{5hv). (3.15) 
Jo 

Therefore, substituting (|3.15p into (j3.14p . one has 

aijdj{5h^)djvdx = 0. (3.16) 

Take v = if5h^ in (l3T6]l . where ?y G (^|rJ , = 1 in ^h, \Dr]\ < — . Then, 



= / aijdj{Shp)dj{7]'^6h(p)dx 

rf'aijdj{5hp)dj{5hp)dx + 2 / aijdj{6hp)r]djr]5h'pdx. 
B2R J B2R 

It follows from Holder inequality and ()3.17p that 

t]'^aijdj{6h(p)dj{6h(p)dx = -2 / aijdj{6h(p)ridjr]Shipdx 
B2R J B2R 

1 1 

2 / /■ „ \ 2 



<2 / r] aijdi{5hp)dj{5hip)dx \ / aijdir]djr]{6h(p) dx 

\ JB2R J \ JB2R 

namely, 

r]'^aijdi{6h(p)dj{5hip)dx < 4 / aijdir]djr]{5h^fdx. 

B2R. J B2R 



(3.17) 
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Consequently, by the strong ellipticity of Ojj , one gets 

|2. 



A / \V{6hf)\ dx < aijdi{6h(p)dj{dh'p)dx 
< 4 / aijdir]djr]{6h^fdx 

JB2R 



and then 



/ \V{5h^)\^dx<C^ I \V^\^dx, yh<R. 



(3.18) 



Therefore, according to (|3.18|) and regularity of minimizer 99, we can conclude that V^(/9 G 
L^{Br). □ 

Next, we derive the boundary estimate of the minimizer ip. 
Lemma 3.5. (Boundary Estimate) For any xq G dVl]^i2, 

vV € (bso (xo) nnA ■ (3.19) 



Proof. Set C4o,5o = ^5o(a^o) n 0^, and 

dxi dy 



Txo ■■ Uxo,5o ^ B+-x^y, y = T^oix), (7ij{y) = -^{v), J{y) = 



u 



For simplification, we write Ux^^Sq and i?^ as U and B^ respectively in the remaining of the 
proof. 

Then for any i/j € Hq{U), ip = il; o , 

6(|V(^p)V(^ • VV'fix = / Q{\aai3di3'f\'^)aijdj(paiidiipJdy, 

Jb+ 

where (p = Lp o . Taking tp as the k-th. difference quotient 

6-hi> =^ ^ (i'iy) - i'iy - hek)j for A; = l, 2, • • • , n - 1, 

we may get from the property and the "integrate by parts" formula for difference quotient that 
for suitable small h > 

= j^^ 5h (^Q{\crai3di3(p\^)aijaiiJdj^pj Biiidy 
4 [Q{\(Tai3df5(f\^)aijdj(p) auditpJdy 



Set 



B+ 

+ j ^{Q{\aaj3dj3(p\^)aijdj(p^ 5h {cruJ) di'ipdy 
= 1 + 11. 

Aj = C Aij{t)dt, Aij{t) = Q{\q{t)\^)6i, + 2Q'{\q{t)\^)qi{t)q,{t), 
Jo 

q{t) = tq\y) + (1 - t)q{y) = {qi{t) , q2{t) , ■■■ , (?„(*)), 
9^(y) = crai3{y + hek)di3ip{y + hsk), q{y) = <yap{y)dpp{y). 
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Now, the term / can be rewritten as 

4 (Qi\(7ai3di3<f\'^)crijdjLpj audiijjJdy 

j ^ Aij^ [q^ - qj^ audii^Jdy 

Aijajsds{6h'f)c7iidii>Jdy+f Aij (j)s'fj Shi(Jjs)criidi^Jdy 



Set 



and 



Then 



B+ JB+ 
h + h. 



tp = ffuh, Uh = Sh'P, 
V G C^{B+), fi^l in B+ = T.,, (u^^ s^ n 17^ ) , |Vr/| < 2 



■0 = -0 o = T^^uh, Uh = Uh o , 
V = VoT,,€C^{UnnL), r/ ee l in C/ nf^L, |V77|<2. 



AijajsdsUhf^iidiipJdy 

B+ 

AijdiipdjUhdx 

u Ju 
hi + h2, 



and 



Aij{ds'f)^6h{(Tjs)(TiidiUhfi'^Jdy + 2 Aij{ds'f)''dh{(7js)(^iiUh'ndif]Jdy 
B+ Jb+ 

= hi + l22- 

Due to the strong elhpticity, 

111 = y Aijri^diUhdjUhdx>\\\viVuh\\\2(i,)- (3.20) 

To estimate the term In, we wih deal with /12, I21, I22 and // first. 
By Holder inequality and the strong ellipticity of Aij, we have 



Il- 



ia 



i_ 1 



2 / AijTjUhdjUhdi-qdx 

<cii(k\\V7^\\l^ - WuhWl^i^u) 

<c4(A||V^||i.(^) 
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and 

1^21 1 = 



Aij{ds(p) 5h{cjs)ciidiUhfi Jdy 

B+ 



< ( / ^AijaudiUhfTjsdsUh'ff'Jdyj ^Aijfi'^6hi(Tjs)ids(p)^dhicrii){di^)^Jdy 



1 



1 



<c/A(A|K|2||V(^iii.(^);' 



1 



<-/n + C^A|i^n|V(^||i.(^)^ 
where K is C^'° norm of the boundary (See assumption ()1.6p ). 

h2 = 2^^ Aij{ds^f6hi(Jjs)(TiiUhvdivJdy < CAK\\\/(p\\l2^u)- 
Next, we estimate //. 

II = j ^ (^e{\aaf3df3ip\'^)aijdj(p^ 5h{(TiiJ)di{'rfuh)dy 

= j ^v"^ (@i\aai3di3'f\'^)crijdjip^ ~5h{cFiiJ)diUhdy 

+ 2 y ^ ?? {Q{\aapdp(f>\^)aijdj(^ 6hi(JaJ)uhdifjdy 

= Ih+Il2. 

Then direct computations yield that 

IIi = j ^ fj^ {Q{\cFai3dj3(p\^)aijdj(^^ 5h{(TiiJ)diUhdy 

<CAK [ f\V^\\Vuh\dy 
Jb+ 



< CAK [ i]'^\V^\\Vuh\dx 
Ju 

<j\\riVu^\\h^^^ + jA'K^V^\\l^^^, 



and 



Ih < 2AK / \V^\\yfi\\uh\dy < CAKWV^Wi^^^y 



Therefore, noticing that 

Q = h+l2 + Ih+ Ih = hi + h2 + hi + h2 + Ih + Ih, 
and applying the estimates as above, we get 

2(C/) -j\\'n^'^h\\\2^jj) + yVi II V f^||i2([;) 



In < CA{K^ + l)||V^||i.(j,. + -hVuH\\l2,u. + -A^K^\\V^\\l,, 



Then, combining with (j3.20p . we obtain the gradient estimates for the fcth difference quotient 
Uh (A; = 1,2, - • • ,n - 1), 

Furthermore, the following derivatives estimates hold, 

n-l 

E < C{K^ + mV^Wl^^uy (3.21) 

fc=l 
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For the D^n'P^ the potential equation, and the estimates for DIj, 1 < k < n — 1, I < j < n, 

%(^B+) - ^(^^ + 1)II'^V'IIl2{;7)- (3.22) 

Combining the estimates ([320, and the F^-estimate ([MD yields (IXT9]l . □ 

Proof of Pro p osi tion \3.3\ It follows from Lemma [331 Lemma 13.51 and a finite cover argument. 
□ 

3.4. Local average estimate. Set 

^xo,r = ^ln{x = {x', Xn) ■ \xn - Xo,n\ < r}, where Xo = (Xg, Xo,n) G ^■ 

Proposition 3.6. (Local CLvevcigc cstiniate). For any xq (£ with \xQ^n\ ^ 77-^' one has 



I a, 

where C does not depend on xo,L 



— I \Vip\^dx < Cml, (3.23) 



Proof. For any < a - 1< a < 6 < 6 + 1< ^, define r] G C°°(Ol), 0<r]<l, |V??| < 2 by 



r]{x) = < 

For any constants ki,k2, set 
ifix) -- 



0, x„ < a — 1, 

1, a < Xn < b, 

0, Xn>b+l. 



ip{x) - ki, Xn<a, 

(p{x) - ki 1 -{xn-a), a<Xn<b, 

— a 



y{x)-k2, Xn>b. 
Then r]'^0 € H^{Q.l) and {jf'p) |a;„=±L = 0. Therefore rj^ip € and 

/ e(|V(/jp)Vv9- V(7?2^)dx = 0. 



Thus, 



/ 7?^e(|V(^p)V<^ • Vc^dx = -2 / r]e{\Vip\^)V^-VT]^dx, 

k2 — ki ^ ^ 

where V99 = Xa,b{x)en, e„ = (0, ...,0, 1), 

— a 

^a,b = {x = ixi,X2, - ■ ■ , Xn) G < Xn < b} 

and Xa,bix) is the characteristic function of ^a,b- Then, 

/ v'ei\^^\')\^^\'dx + / r?2G(|V^p)|^ f-T^) 

2\ 



r?e(|Vv3r)Vv3 • Vrjifdx. 

-1.6+1 

Since ?y = 1 on and / Q(\Vip\'^)— — dx' = mo, 

Jsx^ dxn 

I Tfe{\V^\^)\V^\^dx = -2 I r]e{\Vip\^)Vip-Vr]^dx + {k2-ki)mo. 
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Consequently, 



< 2 

< 4A 

< 4A 

+ 



+ 



6,6+1 



+ \k2 - ki\mo 



+ 



+ 1^2 - h\mo 
\ip — ki pdx 



\Vip\^dx 



6,6+1 



\(f — k2^dx 



6,6+1 



+ 1^2 - ki\mQ. 



Set 



ki 



(fdx, k2 = + (fdx. 
It follows the uniform Poincare Inequality that 
f \ip-ki\'^dx<C 1 \V^\^dx, 

•^^a — l,a "^^a — l,a 

where C does not depend on a, 6. 
Therefore, substituting ([3:25]) into (l3:2i]) yields 



|(^ — k2'^dx < 



6,6+1 



C I \V^\^dx, 

f^6,6+l 



A / \V^Ydx < CA 
We now claim that 



f^a-l,6+l\^a,6 



\V(p\'^dx + \k2 — ki\mo. 



f^a-1,6+1 



\S/ip\dx. 



\k2-ki\ <C 

Assuming (|3.26p for a moment, one gets 

[ \Vip\^dx<C^ [ \Vip\'^dx + ^ 

J^a,b -'^^a-l,6+l\f^a,6 



mo 



f^ti-1,6+1 



\\/Lp\dx. 



Set 5max = max \Sx„\, -S'mm = inin \Sx,^\- On another hand, 



mo 



f^a-1,6+1 



1_ 

\Vif\dx < moSmax 



f^a-1,6+1 



\\/ip\^dx] {b-a + 2) 



<e [ \v^\^dx + ^{b-a + 2)ml 

Combining (13:271) and (l3:28]l leads to 

/ \V^\Hx <Cj [ |V(^| 

•^^a,b ^a~l,b+l\^a,b 



,2j C 

' dx + — e 



f^a-1,6+1 



iVfl'^dx 



C 

+ — (6-a + 2)mg. 



Taking = — yields 



I [ |V(/.|2dx < (Cj + 1) [ \Vip\^dx + ^{b-a + 2)ml 

^ J^a,b \ ^ J^a-l,b+l\^a.b 
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Therefore, one has 

[ \V^\^dx<(c^ + l] [ \\/f\'^dx + ^{b-a + 2)ml, 



16. 



4- 

yA 



/ \Vv\^dx<(c^ + l] [ \V^\^dx + ^ib-a + 2)ml. 



A 
A 



Set = ^ ' . Then < 6lo < 1 and 
C4 + 2 



/ \Vip\^dx<eo[ \Vip\^dx + ^ (6 - Q + 2)mg. 



(3.29) 



Set 



It follows from ([3:29]) that 

6-a+2 C b-a+2 2 

Aa,b < 00-^3^ A-i,.+i + A(CA + 2A) b-a "'O' 

1 + ^0 

Taking = — ^ — and a positive constant A;(0o) ^ 2 such that, if 6 — a > k(0Q), one has 

— a — a 



Then, 



Aa,b < O'oAa-i^b+i + X(cXT2A)'^° yb- a> k{9o). 



A(CA + 2A)(1 
Applying (|3.6p . one has 



/ nAT |^a-Af,b+Af I ^ 2 ^'m, 



6-o + 2iV A(C7A + 2A)(l-0^) 

— '-^I'^O/ '-'maxi'i-Q ~r 



A(C7A + 2A)(1 -0^)' 

2 2 

Therefore, for any — L < a < b < —L and b — a > k^Oo), letting N ^ oo yields 

3 3 



Aa,b < Cml, 



where C does not depend on L. 
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Then, for any xq, |3;o,„| < ^, 



< 



1 



^xo,l| 2A;(6'o)5'max 



< 



2fc(^o)5'max^ 2 
rp- i Gmn 

I 2:0,1 1 

2fe(6'o)S'max 2 
< C IUq, 



which yields ([323]). 

Now, it remains to prove the claim (|3.26|) . For any a G R, we define 



a,; 



^a + i — l,a^ 



i+4 



i-\-i — l,a-\-i 

By the uniform Poincare inequality (j2.3p . one has 



/ 



i = 1, 2, • • • , n. 



^ a — 1 , a 



93 — Qojc?x<C / iVi^ldx < C / \Vip\dx, 



^ a — 1 , a 



S^a-l,a + l 



/ 


(f — ai 


dx < 


(/9 — ai 




2 


f^a-l,a + l 


2 



dx<C ! \Vip\dx. 



Then, it follows (IHTBD]) and (I53T]) that 



Consequently, 



/ 


Q! 1 — CKo 


dx<C [ 




2 


Jn 



S^a-l,a + l 



I V(/9|dj;. 



Q;i — Qo 
2 



< 



C7 



a-l,a| JQa-l,a 



\Vip\dx < 



c 



+1 



\S/ip\dx. 



mill J^a-i „+i 



Similarly, 



Hence, 



«! — ai 

2 



< 



C 



5, 



mill Jfi 



\V<p\dx. 













ai — ai 


+ 


a 1 — ao 


< " /■ 




2 




2 


'S'min Ju 



\Vip\dx. 



In a similar way, one gets 



C f C f 

\a2-ai\<- — / \Vip\dx, \an - an-i\ < — / 

'-'min JQ,a,a + 2 '-'min JQ,a+n-2,a + n 



\VLp\dx 



Therefore, by induction, it holds that 



Wn - aol < 



C 



min Jna-i,a4 



\V(p\dx, 



which proves the claim ()3.26p . 
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3.5. C^'" regularity of the weak solution. 
Lemma 3.7. (Gradient estimate). It holds that 



(3.32) 



where C does not depend on L. 



Proof. The proof is based on Moser's iteration technique. 

Step 1. Interior estimate: It fohows from the definition of weak solutions that for any 
B2R c 

/ G(|V(^p)V<^ • V^dx = 0, VV- G Co~(i?2R). 

'JB2R 

Regarding dgtp as a test function, s = 1, 2, • • • , n, one gets 

0=/ e{\V^\'^)Vip-V{dsi^)dx = - [ ds{Qi\V(p\'^)Vip) -V^pdx 



B2 



Bo 



B2R 



B2R 



Q{\S/ip\^)5ij + 2e'{\Vip\'')d^ipdjip]d,{dsip)dji;dx 
aijdiWsdjipdx, 



where 



= ei\Vip\^)5ij + 2e'{\V^\^)d,ipd,ip G L^iB2R), Ws = ds^ G L'{B2r). 



Therefore 
Taking 

in ()3.33p . one has 

= 



/ OijdiWsdjil^dx = 0, e HI{B2r). 

JB2R 

^ = 7]^wP-\ r]eC^{B2R), r? = l in Br, p > 2. 

aijdiWgdj {rj^w^s'^) '^^ 
B2R 

{p — 1) / rf'aijdiWsdjWsW^g~'^dx + 2 r]aijdiWsdjr]w^~^dx. 



(3.33) 



B2R 



B2R 



Therefore 



{p — 1) / rf OijdiW sdjW sW^g "^dx < 2 
J B2R, 



B2 



< 2| / rfaijdiWsdjWsW^g "^dx 

\''B2R / \JB2R 

1 /* 2 

2 Jb2r P-1 



rjaijdiWsdjTjw^ ^dx 

ttijdiijdjrjw^dx 



aijdirjdjrjw^dx, 



16. 



' B2R 

Due to (13.411 ■ we have 



rj'^aijdiWsdjWsW^ "^dx < 



{p-iy 

4A 



B2 



ri'\Vws\^wP-^dx< . 
B2R (P - 1)^-^ Jb2 



B2R 

aijdiTjdjrjuPgdx. 
\Vr]\^wldx. 



Since 



2|v7 |2 n-2 * 



2 1 2 ' 



V I r^wi — Ws Vr] 



2 2 
> - 



V ( r]w^ 



2 4 



(3.34) 
(3.35) 
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Combining (f33i]) with ([3:35]) yields that 



V ( 7]WS 



4 f 2 - 4A 
dx < ^ |Vr/| w'^dx+ , 



\V7]\'^wPdx, 



namely, 



Bo 



dx <2 



(p-l)2A 



+ 1 



Bo 



\Vr]\'^wPdx. 



Then the Sobolev's inequality implies that 



B2R 



dx 



< C 

< C 



B2R 



V ( rjw. 



dx 



(p-l)2A 



+ 1 



(3.36) 



\S/i]\^wPdx. 



B2R 



for s = 1,2, 
Set 



, n. 



Pk=P 



n 



n-2 



Rk = R{l + 



2k / ' 



C 



Rk — Rk+1 R 

Note that, {pk} is a strictly increasing sequence and tends to infinity as k ^ +oo, and {Rk} is 
strictly decreasing sequence and tends to R as k goes to infinity. The following is the standard 
Moser's iteration process. 

Taking p = pk, rj = rjk in ()3.36p yields that 



w7^'dx < C 



piA 



< C 



{Pk - l?\ 
K2^ 



+ 1 



Br^ 



XR 



and so. 



Let Mfc = 



c 



A2^ 

Tr 



and Dk 



C- 



AR 



. Then by induction 



Mfc+i < DfcMfc < • • • < L'fcDfc-i • • • DoMo = Mo [] Dj. 

j=0 



(3.37) 



Due to 



n«.=n 

i=o j=i 



c 



A2^ 

Jr 



c 



A_ 
XR 



2--j=up, = c C 



A \ 2p 



so taking /c — )> oo in ()3.37p yields 



A \ 2p 



supw;^ < C|^C— j \\ws\\lp{B2r)^ for any s = 1,2, 



,n. 



(3.38) 



Step 2. Boundary estimate: For any xq G 90^/2) according to Lemma 12.51 there exists a 
neighbourhood C/^o of in M" and an invertible C^'" map 

Txo : C4o 1^ ^L/2 ^ Bj^,-x^y 
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satisfying ()2.4p . where Bs^^ is independent of xq. Define 



Then ([23]) implies 

o-jj(2/)o"m(y) = 0, on ?/„ = 0, j = 1,2, • • • ,n - 1. 

For any < R < 6o, 

1 eaijdj^aiiditl^Jdy = 0, Vtl^ e hUb+] . (3.39) 

Taking (s = 1, 2, • • • , n) 

i^ = ds (r?2 (a,^)""') , fi€C^ (bI) and 7? ^ 1 in B+ 
in (|3.39p and integrating by parts show that 



- Ids (^QcTijdjifxJiiJ^ di {rf{ds'^Y~^^ dy 

/ [QcTijdjipaii) di (r?^(4(^)*'~M J6sndy' 

JdB+n{y„=o} ^ ^ ^ ^ 

- / 4 (Qaijdjipauj) di {7f{ds(pY'^\ dy, 
J si 



(3.40) 



where the boundary terms vanish according to (|2.4l -2) and (|3.13p . 
Detailed calculations show that 

ds{Qcrijdj(faiiJ) = Qaijdj{ds^)cFiiJ + 2Q' aai3dp(pds{cFa'^d^^)aijdjLpaii J 
+ QdjipdsioijauJ) 
= @aijdj{ds(f)aiiJ + 2Q'aai3di3(paa^d.y{ds<f)aijdjifaiiJ 

+ 2Q'aai3di3(pds{cra-y)d.ytfaijdj(paiiJ + QdjifdsiaijauJ) 

+ 2Q' aai3dpLpds{(Ta^)d^i^aijdji^aiiJ + Qdj(fds{aijaiiJ) 
= i&Sia + 2Q'aai3di3(paijdj(f)aa'yd^ids'f')(7uJ 

+ 2Q'aai3di3(fds{cra^)d.y(paijdj(faiiJ + Qdj(pds{aijaiiJ) 
= {Q5ij + 2Q' aiada'f>cFjiidp(p)ai^d^{ds(p)cFjiJ 

+ 2Q'aa/3d/3ifds{aa-y)d^'f(Tijdj(paiiJ + Qdj(pds{(TijaiiJ). 

Set 

Ai-y = [Q5ij + 2Q'aiada'f(7j(sd(3(p)ai^ajiJ, 

Bis = 2@'aapdi3(pds{aa'y)d^(faijdjtfaiiJ + Qdj(pds{cFijaiiJ). 

Then 

! Ai^~d^{dsm{r}\ds^Y~^) + Bisdi{f{ds^y-^)dy = 0. (3.41) 
Jb+ 
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Denoting ds(p by Ws, we have 



which can be rewritten as, 



Note that, 



Jb+ Jb+ 

+ {p-l) I Bi,fwl"^diwsdy + 2 / BisfjdifjvjP-^dy 

^ R ^ R 

= I^ + l2 + l3+ h. 



\h\ < 2 



fiAi^difjd^Wswl ^dy 



< 



p — 1 



Bt 



fAi^diWsd^Wswl "^dy-l / Ai^difjd^fiwPdy, 

P - i Jb+ 



and 



rfSisdiWsw'^s "^dy 



Bl 



<{p-iy- I fj^\Vw,\-'wP-^dy + ^ I f\Bis?wl-^dy, 



A 

'4 



I/4I < 2 



vBisdifjujP ^dy 



Bt 



A .Ib+ 



<2 f,\Vfi\\Bis\wP~'dy. 

B+ 



Therefore, substituting I^Mh . and I^M\i into <^M\i yields that 

-^1 _ {P 1) /" fj^Ai^diWsd^Wswl~'^dy 



< I Ai^BifiB^fjwldy +{p-l)j [ f\VwsfwP-^dy 

P-^ Jb+ 4 Jb+ 

+ ^/ fl'\Bis\'wr'dy + 2 [ fj\Vf,\\Bis\wP-'dy. 
A Jb+ Jb+ 



Then the uniform ehipticity yields 



Bl 



f]\Vfl\\Bis\wl~^dy. 



Note that p.3p implies that 



\Bis\ < C(l + A)|u;| < CK\w\, 
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where w = (t^i, ^2, • • • , Wn) = {di(p, d^^^ • • • , i9n<^) and = \'Ws\- Then, 



Bz A(p - ly Jb+ 



\Vf,\^wPdy + 



CA2 



~2| ~ |2 ~p-2j 

fj \w\ w'g ay 



which imphes 



[ f V (wi 

Jb+ ^ 



dy< 



+ 



CA 



A(P-I) JBt 



i]\Vfi\\w\w^ ^dy 



A(P-1)2 JBt 



\Vfi\^wldy + ^p^ / f\w\^wP-'^dy 



Bl 



+ 



It fohows from ([OT]) that 



< 



A(P-1)2 



A(P-I) 75+ 



+ l) / \Vf,\''wldy + ^p' [ fi^\w\^wl-^dy 
J Jb+ a Jb+ 



for s = 1, 2, • • • , n. 
Let 



1 



Pk=P 



n 



% G Co"^(Si?,), % = 1 in 5/?,+, and |Vr?fc| < 



/c = 0,1,2,..., 

2 2^^'+^ 



Rk-Rk+i {l-e)6o' 
Taking p = pi^, fj = fj^ m ()3.48p and using Sobolev embedding Theorem, we obtain 



w^'^'dy 



Bt 



< I I [ VkWs 



2n n — 2 

Pk\ n-2 
2 



dy 



< C 



Bl 



V fikw. 



Pk 



dy 



+ 



KPk - 1) Jb+ 



flk\^fik\\w\wl'' ^dy 



-^Wpk-iy^' 



4fc+i 



CA f 
w^'' dy + —T^Pk^ I \w\^wl''~'^dy 



(3.47) 



(3.48) 
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Therefore, 



\LPk- 



,(^+ \< (Ak [ wP^dy + Bk f \w\^wP^-^dy + Ck [ 



iwlw 



Pk- 



with 



and 



Ck 



■Pk'i 



\{Pk-i) {i-e)8o 



Note that 



/ w^'' ^\w\dy < 
Jb+ 



Bl 



w^'' dy 



Pk-2 
Pk 



Pk 



Bt 



B7 



\w\'P'' dy 



\w\'P^ dy 



2 

Pk 



1 



Then 



||Pfc-2 
LPk ( B 



\W\ 



'LPk B 



+ Ck\\Ws 



\w\ 



< \\w. 



Pk 



'LPk{B^) 



)+Bk\\\w\\\lvk(B 



Rk> 



+ Ck\\Ws 



\LPk(B+^)\\M\\LP^B+^) 



1 

Pk 



Hence, 



Ak + Bk + Ck 



1 

pfe 



2 

I Pk 



LPk B 



which imphes that 

n 

E 



< 



1 

Pfe 



< 



Afc + Bfc + Ck 
{Ak + Bk + Ck)n^ 

n 



\W\ 



2 

I Pfc 

LPfe ( Bt 
2 

I Pfc 

LPfe ( B 



E 

■^fc/ s=l 



Pfc- 



Pfc 



LPfc B 



«fcy 



" Pfc -2 Pfc 

\\w II '^^'"^ 



Rk 



n 

E 

pfc-2 
Pk 



'LPk Bt 



2 

Pfc 



Therefore, 



< 



s=l 



'^k + lj 



{Ak + Bfc + CkW 



E 



lLPfe(i?+ 
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Define 



Then 



■w. 



'LPk I B 



It is clear tliat pk = p 



s=l 



n 



Hi 



n-2 
Ak + Bk + Ck < C 



Mk+i < DkMk. 
< p^ , for n > 3, so 



2Kpl 



#+1 CA2 2 CKpl 2'=+! 

+ 1 I — ^sn H T^Pk + 



Krt -iV- ; (1 - fi?il A2 ' \{pt - 1) (1 - 9)* 



4^ 



A ){l-eY5l A2 



A"(l-e)<5o 



where T = C 
Then, 

Note that 
One has 

Letting k ^ oo shows that 



<T ■ l6^ 

1 A A2 2 A P 



(1-0)252 A ^ A2^ + A {I -6)8, 



C does not depend on k. 



Mk+i < DkMk < Dk ■ Dk^i ■■■Do- Mo 



El n \ - ? 

— = — , and > — < 00. 
Pi 2p ^ Pi 

1=0 2=0 



Mk+i < CT^Mo, V/c > 0. 



sup ItZ-sI < CT^ 

S = l ^0^0 



'LP B 



(3.49) 



Combining the interior estimate p.38p with the boundary estimate ()3.49p yields the desired 
gradient estimate p.32p . □ 

Remark 3.2. It has been assumed that Wg > and Ws is bounded in the above proof. The 
boundness assumption could be eliminated by a standard technique (see chapter 8 of [M]). If 
Ws is not positive, we can repeat the proof for wf and respectively. 

Remark 3.3. In the case that n = 2, choosing pi = 00, one can obtain the estimate similarly to 

Lemma 3.8. (Holder estimate of gradient.) ip G C^'° o,nd 

l|V</'llc70."(n^/,) ^C-mo, (3.50) 

where C does not depend on L. 

Proof. Step 1. Interior Estimate. For any B2R C ^, Wg = ds(p (s = 1, 2, • • • , n) is a weak 
solution to 

di{aijdjWs) = 0, 

in the sense of (|3.33p . where aij = p{\V (p'\^)5ij + 2p' (\'\/ (p'^)di(pdj(p. Then, the desired interior 
Holder estimate for Ws is just the standard interior Holder estimate for the weak solutions to 
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second order elliptic equation with bounded coefficients. 

Step 2. Boundary Estimate. Similar to (|3.4ip . one has for any s = 1, 2, • ■ • , n — 1, 



Ai^d^wsditp + BMdy = 0, e ^^o(^A) 



(3.51) 



where B^, Aiy and Bis ai'e same as in ()3.4ip . 

By an even symmetrizing procedure, Ws, Ai^ and Big denote the even extensions of Wg-, Ai^ 
and Bis, respectively. Then Ws satisfies 



/ Ai^d^WsBd + BMdy = 0, HI{Br), 1 < s < n - 1. 

Jbr 

Since ||u;||loo < Cmo, 

11^75 -BjsIIloo < C, \\Ai^,Bis\\l°° < C. 
Therefore, for 1 < s < n — 1, the standard interior De Giorgi estimate gives 



(3.52) 



+ \\\Bis\\T,qiB+\ ) < CiUQ, q> n. 



Now, we estimate Wn- For any yo G — | 6^' 2 }' ^^^^'^^ ^ ~ v'^i^s — Ws) in (|3.52|) . 

?7 G C^(B2r-(2/o)), r/ = 1 in Br{yo), |Vry| < - and Wg = f Wgdy, one has 

^ JB2r{yo) 



Ai^d^Wsdi{fj {lis - Ws)) + Bisdi{fj {ws - Ws))dy = 0. 



Therefore, 



< 



< 



Bl 



Ai^diWsd^Wsifdy 



BisdiWsffdy 


+ 








Jb 









Ws)dy 


+ 


2 / Bis0ifj{ws - Ws)dy 






Jb+ 



A 



^5 
4A 



Bis\^fi^dy + - / \Vwsrfi'dy + 
4 JBt 



4A 



Bl 



Ai-ydiWsd^Wsrfdy 



^ Jb+ 
As a consequence, 

ff\Vws\^dy < C 



+ -T- I Ai^difid^7){ws - WsYdy + I \Bis\^Tfdy + 



Bl 



{ws - Wsfdy. 



Bl 



Bt,r\B2r{yo) 



\Vfi\^\ws - Ws\^dy + 



B+nB2r(yo) 



f\Bis\^dy]. 



Noting that \\ws 



R/2) 



< CniQ and \Bis\ < Citiq, one has that, for any 1 < s < n — 1, 



B^nBrivo) 



|Vu),p < C 



B^nB2r{yo) 



\Vfi\'^\ws - Ws\^dy + 



f\Bis\''dy 



B^nB2r{yo) 



< Cml (r"-2+2" + r") 



(3.53) 
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According to the equation of c^, one has 



JB+nBr(_yo) \ ^^iJB+nBriyo) J B+nBr{yo) 



Then, due to ()3.53p and (|3.54p . one has by Theorem 12.11 that 



(3.54) 



Now, the Holder estimate of Vip follows from Step 1 and Step 2. □ 

3.6. Proof of the existence of subsonic flows. Proof of the statement (i) of Theorem 

1. For any fixed suitably large L, according to previous subsections, one can get a function 
such that [ipi{x) — (/3l(0)) G Hl is a weak solution to problem [2j Set (pl{x) = ^l{x) — 
(Pl{0). Moreover, ipi G (^^'"(n^/a) and 

||V(^L||^o,a(ni/2) < Cmo. 

For any fixed K » 1, if L > 2K, 

ll<^L|lci."(nK) - 
where C does not depend on L, and (pL satisfies 

/ e(|V(^L|2) V(^i • vvdx = o, y^jec^inK). 

Since ipi ^ Hi n C^'"{ftK) satisfies the equation (|3.10p . one can check easily that 

/ G (|V(^Lp) — dx' = mo, for any xq G ^k- 

By a standard diagonal argument, there exists a 99 G C^''^{^1) and a subsequence (fL„ such that 
for any K, 

lim ll'^Ln - '/5|lci."(!^x) = 0- 
n— >cc ^ ' 

Therefore, one has 



and 

It is clear that 
and 



/ e(|V99p)V(^ • vV'dx = 0, vV'gc^(o), 

Jn 

j 0(|V(/9p)^ — dx' = mo, for any xq E il. 

J SxQ ^ 

<^GCi'"(o)n^^L(^^) 



l|V<^IIco.-(f7) < CniQ. 

Similar to the previous subsections, one can prove that (p G Hf^^{^) and 99 is a strong solution 
to 

' {Q{\Vip\^)5ij + 2Q'{\Vp\^)d,pdjp) d^if = in f], 

< ^ (3.55) 

— ■ = 0, on oil. 

By the standard regularity theory for second order elliptic equations, one gets that ip G Cf^^{Vl) 
is a solution to ()3.55p with the property 

l|V</?||ci.-(n) < CrriQ. 
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Choose mo small enough such that Cttiq < 1 — 26q. Then the subsonic truncation automatically 
disappears, so ip & ^loci^) ^ smooth solution to the original Problem [TJ This proves the first 
part of Theorem II. li 

Remark 3.4. In fact, we can conclude that (p € C°°(r2) by the standard bootstrap argument. 



4. Uniqueness of the global subsonic flow 



Theorem 4.1. (Uniqueness) Suppose that satisfies the assumptions il.6\) . and (pk {k = 1,2) 
are uniformly subsonic solutions to the following problem 



diu (/9(|V(^fc|^) Vi^fc) = 0, in Q 



0, 



. dfi 

associated with the same incoming mass flux rtiQ. Then 

Vipi = S/ip2, in ^■ 

Proof. Set (f = ipi — (p2- Then ip satisfies 

di{Aijdjp) = 0, in 



on 90, 



dp 
dfi 



0, 



on 90, 



where 



^ij = [ Pi(f)^ij + '^p'{Q^){sdjPi + (1 - s)djP2){sdiipi + (1 - s)diip2)ds, 
Jo 

= \sVpi + (1 - s)Vp2\^. 
Moreover, there exist two positive constants A < A, such that for any vector ^ G 

Let r/(x) = 7]{xn) be a function satisfying 

Vi^n) = 1 for \xn\ < L; rj{xn) = for > L + 1, and |?7'(x„)j < 2, 



h = {x = 


(x , x^ 


) G $7 a 


< Xn < b} and for L > 




p{x) 


- V'Z' 


X G n {Xn < -L}, 


(p(x) = < 


p{x) 


-^~L- 


^^^^^^{Xn + L), xenn{-L<Xn<L} 




p{x) 




X G n {x„ > L}, 



where 



1 



Note that Vp = Vp 
function of l. 



p{x)dx, p'l 



1 



10 



p{x)dx. 



(4.1) 



(4.2) 



2L 



X-L,L{x)en, Cn = (0, • • • ,0,1), X-l,l{x) is the characteristic 
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Multiplying on the both sides of the first equation in ()4.ip by rf^-, and integrating it over 
r^, one obtains 

/ rf Aijd^^dj^dx + f rj^ (p(|V(/Pi|2)V(/Ji - p{\V ^2\^)V ^2) ■ Sndx 



i]{(f — (p^)Aijdirjdj(pdx — 2 / ri{(p — LpJ^)Aijdi'r]dj(pdx. (4-3) 
The second integral on the left hand side of (j4.3p vanishes. Indeed, 

r]\p{\Vipi\'^)Vipi - /3(|V(^2p)Vv?2) • endx 

n-L,L 

v\t) [ (p(|V(/Jip)Vv9i • en - p(|V¥P2p)V(^2 • 4) (ix'dt = 0, 

-L JSt 

since the two solutions possess the same mass flux mo, here = H {x„ = t} for t € [—L, L]. 
It follows from K3\i and (IMI) that 



< 4A / 1^9 - (/?;f;||V(/?|(ix + 4A / |(^ - (/7;^||V(^|(ix 

< 2A I / |(^ - (/?~p(ix + / |v3 - + / 

Due to the uniform Poincare inequality, ie. 

/ \ip - ifll'^dx <C \Vip\'^dx 

/ \ip-^l\'^dx <C 1 \V^\^dx, 
where C is independent of L, we have 



and 



A 



/ \V^\^dx <C I \V^\^dx. (4.4) 



By the estimate (|3.6|) . one has 

|2 



[ \V^k\''dx<CmU\^-L-i,-L\ + \nL,L+i\), A: = 1,2, 



which implies that 

\Vip\^dx < C, 



I 



-lU ^l,L+1 

where C is independent of L. 

Combining this with (j4.4p shows 

/ \Vip\'^dx ^ and / \Vip\^dx 0, as L 

Taking L — )■ 00 in (|4.4p yields 

V(/3 = in Vt. 



00. 



□ 
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As a direct application of the uniqueness, we can obtain the exphcit form of the subsonic 
solution ip{x) to the Problem [U provided that the nozzle is a cylinder. 

Corollary 4.2. (Cylinder case) Suppose that Q is a cylinder, that is, Q = S x {—oo, +00), S is 
a n — 1 dimensional, simply connected, C^'° domain. Then the unique solution to ProblemUl is 
given by 

^ = qoxn + V'O, 

where ipo is an arbitrary constant, qq is a constant defined by 

p{Qo)Qo = Y^- 



\ u = V(/?= (0,---,go) 
' P = Pilo) 



s \ 



Figure 5. Subsonic flow in cylinder case 



5. Existence of the critical incoming mass flux 



In the Section 3 and Section 4, we have obtained the existence of the uniformly subsonic 
flows associated with suitable small incoming mass flux ttt-q and the uniqueness of the uniformly 
subsonic flow. In the following, it will be shown that there exists a critical mass flux Mc such 
that the flow is always uniformly subsonic, provided that the mass flux mo is less than Mc- 

Theorem 5.1. Suppose the nozzle satisfies the basic assumptions lil.6\) . Then there exists a 
positive constant Mc < 1, which depends only on 0,, such that if < rriQ < Mc, then the 
following problem 

(ira(p(|Vv3p)V(^) = 0, in 

— = 0, on af], 

on 

[ p(|V^|2) ^dS = mo 
V Js ol 

has a unique uniformly subsonic solution ip{x) up to a constant satisfying 

Q{mo) 

Moreover, Q(mo) ranges over [0,1) as tuq varies in [0,Mc 



sup \Vip\ < 1. 



Proof. Choosing a strictly increasing sequence {qn}'^=i satisfying lim g„ = 1. Consider the 

n— >-oo 

following truncated problem 



div {pni\V^\^)Wip) =0, 



dn 



in 
on 



(5.1) 



Pn{\V^?) %dS = m, 
ol 
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where 

p{s), if 0<s<ql 



Pn{s) = \ smooth and decreasing if < s < , ^ 

satisfies Pn{s) + 2sp'^{s) > Xn with some A„ > for all s > 0. Let ipn{-',m) solve the problem 
()5.ip and set 

Q„(m) = sup |V(^„(-;m)|. 

We claim that Qnijn) is a continuous function of m. 
In fact, we take a sequence rrij m, it suffices to prove 

\S/ipn{-;rnj)\ |V¥?„(-;m)|. 

Without loss of generality, we assume that there exists a positive constant M, such that 
supmj < M. 

It follows from Section 3 that the solution ipn{-;mj) to the problem ()5.ip with the mass flux 
rrij satisfies the Holder gradient estimate 

II V¥P„(-;mj) ||ci,-(n)< C(M) (5.2) 

and 

II - (^„(0; jnj) ||c2,<.(Q^)< C(M,L) for any L > 0. (5.3) 

Therefore, by Arzela-Ascoli Lemma and a diagonal argument, there exists a subsequence 93„(-; mj^.)- 
(/?„(0; rrijf^) such that for any L > and < /3 < a 

{ipn{-;rnjj - ipniO;rnjJ) ^ ipn{-) in C^'^(f]L) as mj^ ^ m. 

And 93„(-) solves the boundary value problem (|5.ip and satisfies that 

II VvJn Ilci.-{n)< C'(M). 

On the other hand, it follows from the previous sections that there exists a (/?„(•; m) which solves 
()5.ip . We can conclude that 

Vv3„(-) = S/ipn{-;m) 

by the uniqueness. 

Hence, for any L > 

Wipni-^rnj) ^ \/(pn{-;rn) in C^'^(ilL), as nij m, 

which proves the claim. 

It follows from the claim that, there exists the largest Qn > and the smallest 5„ > such 
that 

Qn-i < Qn{m) < Qn, for any m G {mn,Mn). 
Moreover, clearly Mn+i > M„. Set = lim M„. It follows the definition of M„ that M„ < 

n— >oo 

p{Ql{Mn))Qn{Mn) < 1, hence M, < L 

Then we can conclude that there exists a critical mass flux < 1, for any mo < Mc, there 
is Mn such that M„ > "iO) then 

Q(mo) = Qn("lo) < Q'n < 1. 
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Moreover, for any normalized subsonic speed Q € (0,1), there exists some n, such that 
Q G {0,qn), therefore, there exists a uiq € (0, M„), such that Q(mo) = Qn{mo) = Q by the 
continuity of Qn{fn)- 

This completes the proof of Theorem 15. 1[ □ 



6. Properties of the subsonic flow 

In this section, we consider the asymptotic behavior of the uniformly subsonic flows at the 
far fields under the asymptotic assumption p.7|) . 

Proposition 6.1. Suppose that the nozzle satisfies the asymptotic assumption |j. 7| j. Then the 
subsonic flow constructed before approaches to uniform flows at the far fields, ie. 

Vip = (0, • • • , q±), as Xn ±oo, 

respectively, q± are constants uniquely determined by 

p{q±)q± = j^, 

respectively. 

Proof. Assume that (p{x) is a classical solution of 

div(p(|V(^P)V(^) = 0, in n, 

dip 
dn 

[ p{\Vip\') ^dS = mo, 
V Js ol 



0, on 50, 

on 



satisfying 

II ||ci'-(f7)< CniQ. (6.1) 

Step 1. A Special Case. Suppose that n {x„ > Lq} = C/+ x [Lq, +oo) for some Lq. Define 
a sequence of functions as follows 

ipk{x',Xn) = (p{x',Xn + k)xnk, 

here Qk = {{x' , Xn)\{x' , Xn + k) e^, Xn + k > Lq + 1}. 

For any compact set S C and k sufficiently large, it follows from the gradient estimate 
(O) that 



II ||ci'<»(5x[-fc/2,fc/2])< C, 

where C is independent of k. Set ipk{x) = <pk{x) — (pkiO), for any fixed L > 1, if A; > 2L, we have 

II V'k llc2."(Sx[-L,L])< C, 

with C independent of k. Therefore, by Ascoli-Arzela Lemma and a diagonal procedure, there 
exists a subsequence 0kj i such that for any L 

ipk^ ^ po, in C72'/3(S X [-L,L]) with /3 < a, 
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for any compact set S C 5+. Therefore (po solves the following problem 

div{p{\Vip\'^)Vip) = 0, in E+ = S+ x {-oo,+oo), 

dip 



(6.2) 



„^ 0, in dS-i. X (— oo,+oo), 

on 

I P(|V(/^P) ^d5 = mo, 
Js ol 

ip{x) = 0, on Xn = 0. 

v 

Moreover, 

Vipk = V^u ^ V(/.o in C^^^'iS x [-L, L]) for < /?. 
So, choosing 5" = 5+ and L = 2, we have 

II V(^fc - |lcM(5_^xh2,2])^ as A; ^ +00. 

By the definition of ipk and Corollary 14.21 it follows that 

Vip — > V(^o = (0, • • • , q+) as — >■ +00. 

This completes the proof of Proposition 16.11 in this special case. 

Step 2. General Case. Suppose now that the nozzle satisfies ()1.7p . we can also define a 
sequence of functions as 

ipk{x',Xn) = (p{x',Xn + k)xnk, 
here = {{x' , Xn)|(x', Xn + k) & Q, Xn + k > 1}. Then similar to the Step 1, we can show that 

V<^fc^. ^ Vipo in C^'^iS X (6.3) 

for any compact set S C S+ and any fixed L, here 5 may not reach the boundary dS+, and ipo 
is still the solution of boundary value problem ()6.2p . 

In particular, c/^q satisfies the no-flow boundary condition on the nozzle wall. Indeed, for any 
given point (y', ?/„) G dS+ x (— 00, +0x0), n = (ni, 0) is the outer normal direction of the cylinder 
S+ X (—00, +00) at {y',yn)- For any 5 > 0, there exists suitable large Kq > 0, such that 

{y' - 5ni,yn + k) e S X {x„ = yn + k} for k > Kq, 

where 5" is a compact set of i7 n {x„ = yn + k}. 

There exists a sequence of n — 1 dimensional vectors {z^}^]^, such that (y' — 5ni + z^, + 
k) G (90, and |ifc| = dist{{y' — 6ni,yn + k),dQ). is the out normal of the domain at 
{y' - 5ni + Zk,yn + k). Obviously, 

lim \zk \ — >■ 0, and lim n^. n, 

fc— >+oo fc— )-+oo 

due to the assumption ()1.7p on the nozzle at the far fields. 
Therefore 

Vipo{y' - 6ni,yn) ■ ft = lim Vip{y' - 5ni,yn + k) ■ n 

fc— )-+oo 

= lim (V(/9(?/' - (5ni, y.„ + k) - Vip{y' - 5ni + Zk,yn + k)) ■ n 

fc— >-+oo 

+ lim Vifiy' - dfii + 4, yn + k) ■ {ft- Uk) 

fc-^+00 



0. 



As a consequence, -rr^iu' ^Vn) = 0. 
an 
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Set 

n\s) = {x' enr\{xn = k}\ dist{x', an) < 6} , no{6) = {Qr\ {x„ = k}) \n'{6), 

and 

Bi4y'^ = {x' eW'-'\\x' -y',\<26}, y', e On n {x^ = k}, i = l,---,N, 

such that 

N 

n'{6) C U B,4yi), for any ,5 > 0. 

i=l 

For any fixed 6 > 0, there exists a sufficiently large Kq > 0, such that ^q{S) C 5 for A; > Kq, S 
is a compact set in 5+. 

Choosing L = 2, one has from (|6.3p that 

II V(^- Vv^o ||c^(f^,(5)^[fc^;-,+2])^0 as /c^+oo, for /i</3. (6.4) 

Near the boundary dO, fl {xn = k}, ip possesses the following estimates 

II \\cv{B+^{y'^)x{k,k+2))^ 

for ry > 0, with Bf^^{y[) = Bi^siUi) n (0 n {x„ = /c}). Hence, for any x',y' G -Bj^^lyO; one has 

\V(p{x',k) -V^{y',k)\ < 

Then, for any e > 0, there exists 5 > 0, such that 

\V^{x', k) - V(^(y', ^^)| < for any x' , y' € B+^{y[), (6.5) 

and i = 1,2,--- , iV. 

On the other hand, it follows from ()6.4p that there exists iiT > such that 

iVvJ - Vv?o| < for any x £ fi'o(5) x {K, +oo). (6.6) 

Then, combining that (|6.5p and ()6.6p . one can conclude 

|Vv? - (0, • • • , g'+)| < e, for any x„ > K. 

Similarly, one can get the asymptotic behavior as Xn — oo. This completes the proof of 
Proposition 16.11 □ 
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